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We analyze the universal critical behavior at the chiral critical point in QCD with three degen- 
erate quark masses. We confirm that this critical point lies in the universality class of the three 
dimensional Ising model. The symmetry of the Ising model, which is 7,(2), is not directly realized 
in the QCD Hamiltonian. After making an ansatz for the magnetization- and energy-like oper- 
ators as linear admixtures of the chiral condensate and the gluonic action, we determine several 
non-universal mixing and normalization constants. These parameters determine an unambiguous 
mapping of the critical behavior in QCD to that of the 3d-Ising model. We verify its validity by 
showing that the thus obtained orderparameter scales in accordance with the magnetic equation 
of state of the 3d-Ising model. 
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1. Introduction 

A detailed understanding of the phase structure of strongly interacting matter at non-zero tem- 
perature and density is one of the major tasks for non-perturbative lattice QCD simulations. It has 
been conjectured based on renormalization group arguments §\§ and verified in numerical calcula- 
tions ^ ||] that the QCD phase transition with three degenerate quarks flavors is of first order for 
small masses. This also holds for two light quarks and one heavier strange quark, if the strange 
quark mass is kept below a critical value. The line that separates the first order region from the 
crossover region is a critical line on which the nature of QCD phase transition is of second order 
[Q]. It is an open and interesting question how close the critical line passes the physical point, where 
we know that the QCD transition is a crossover [Q]. In fact, recent lattice calculations suggest that 
the first order region dramatically shrinks with the approach of the continuum limit [|], [7|]. 

In this work we restrict ourselves to three degenerate quark masses. Within this theory the 
critical point - the endpoint of the line of first order transitions - has been determined before, with 
standard staggered fermions on A^ T = 4 lattices ^ ||] and on N r = 6 lattices [0]. Our analysis is also 
based on simulations with standard staggered fermions (and the standard Wilson gauge action) on 
Nt =4 lattices. We repeat the determination of the critical endpoint with larger spatial lattice sizes. 
Moreover we focus on the determination of the correct scaling operators and variables that allow 
for the mapping of the QCD critical behavior onto the corresponding universality class, which is 
here the 3-dimensional Ising model. The construction of the operators and fields is done along the 
line of [Q, ||]. 

For our simulation we use an exact RHMC algorithm which we have implemented to run 
entirely on GPUs. I.e., all calculations that involve spinors or gauge fields are done on the GPU, 
whereas the CPU is only used to control the I/O and program flow. The HMC trajectory is done in 
a mixed-precision manner: the force and evolution calculations are done in single precision, which 
is corrected by a Metropolis step in double precision. The expectation values which we obtain 
for the chiral condensate and the plaquette are shown in Fig. [[]. In oder to analyze the finite size 
scaling behavior we have performed simulations on spatial volumes of 16 3 , 24 3 in some cases also 
32 3 . For each coupling j8 = 6/g we have generated between 5000 and 10000 configurations that 
are separated by 20 RHMC updates. For the bare quark masses, we picked values in the range 
of 0.024 < am < 0.034. From Fig. | one can already see that the transition becomes stronger for 
smaller quark masses. 

2. Z(2) critical behavior 

The dynamics and universal critical behavior of the 3d-Ising model is - as in all generic O(N) 
models - characterized by two relevant scaling variables. Close to the critical point we thus can 
assume that the Hamiltonian takes the effective form 

H eg (t,h) = tE + hM , (2.1) 

where t and h are the relevant scaling variables and E,M the corresponding scaling operators. In 
analogy to a classical spin system we will refer to the fields as reduced temperature and external 
field, and denote the operators as energy and magnetization. 
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Figure 1: The chiral condensate (left) and the plaquette expectation values for different bare quark 
masses as a function of the coupling j3 = 6/g 2 . The values of the bare quark masses are am = 
0.024,0.026,0.028,0.030,0.032,0.034 (from left to right). 

Under a renormalization group transformation the singular part of the free energy scales as 

f s {t,h)=b- 3 f s {y%^h). (2.2) 

The dimensionless scale factor b is arbitrary. From here immediately all known scaling and hyper- 
scaling relations follow. Especially the standard finite size scaling behavior (choosing b = LT = 
N a /N r ) of the susceptibilities is given by 

V- l {(5E) 2 )~L a l v , V~ l {(SM) 2 ) ~ L y / V with 8X=X-(X). (2.3) 

Here L denotes the spatial extent of the system and a, 7, V are the critical indices that characterize 
the universality class. In order to map the QCD critical behavior to the Z(2) critical behavior we 
construct the scaling operators from the operators appearing in the QCD Lagrangian by employing 
the following linear ansatz |§] 

E = Sc + ryy , M = \j/\j/ + sS G . (2.4) 

Here So denotes the Wilson gauge action and yy the chiral condensate. Similarly, we linearize 
the scaling fields in the vicinity of the critical point (j8 c ,m c ) by assuming 

t = t ((p-p c )+A(m-m c )), h = h ((m-m c )+B(P-p c )). (2.5) 

The mixing parameters of the scaling operators and fields: r,s,A,B as well as the normalization 
constants to, ho have to be determined from the QCD simulations. Although there is no apparent 
reason, we assume in the following the directions t and h as orthogonal, which relates A and B by 
A = tan{§) = —B, where (p is the rotation angle between the j3 and the t direction. 



3 



Critical behavior in 3-flavor QCD 



Dominik Smith 




Figure 2: The fourth order Binder cumulant of the chiral condensate (left) and the pseudo-critical coupling 
from the peak position of the disconnected chiral susceptibility (right). 



3. Locating the critical point 

The first step to establish the correct scaling operators and fields in the 3-flavor QCD case 
is the location of the critical point. It has has been shown [Q] that a good method to do so is by 
analyzing Binder cumulants. The Binder cumulant of the order parameter is defined by 

B 4 (M) = ((8M) 4 )/((8M) 2 ) 2 . (3.1) 



From the scaling behavior of the free energy ( J2.2| ) it follows immediately that the Binder cumulant 
is volume independent at the critical point. Moreover, its value at the critical point is universal. Fur- 



thermore, we can conclude form (2.2) that close to the critical point it is not mandatory to analyze 
the Binder cumulant of the correct order parameter. As long as the operator of consideration has 
non-vanishing overlap with M, the strongest singularity corresponding to Ul v will dominate and 
the correct universal value of B4 is obtained. This means that in our analysis of Binder cumulants 
we can choose the mixing parameter s arbitrarily 1 . In Fig. || (left) we plot the Binder cumulant 
of the chiral condensate (s = 0) along the pseudo-critical line. I.e. for each mass we evaluate B4 
at the pseudo-critical coupling j8 pf (m) which we have determined from the peak position of the 
chiral susceptibility. The two volumes V = 16 3 and 24 3 clearly suggest an intersection point close 
to the universal value of the 3d-Ising model, which is B^ mg = 1.604. We can thus confirm that 
the critical point lies in the Z(2) universality class. Furthermore, we obtain for the intersection 
point of the 16 3 -lattice with the universal Ising value a mass value of am c = 0.0268(6), which is in 
good agreement with the previously obtained value from RHMC simulations am c = 0.0263(3)[|]]. 
Combining the results of the 24 3 , and 16 3 -lattices into a common fit that restricts the intersection 
point of the two volumes to Ising universal value, we obtain a critical mass of am c = 0.0273(4). 



'As long as we have not 5= l/r, where the leading order singularity exactly cancels. 
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This fit is shown as straight lines in Fig. || (left). For both of the above described fits we restrict 
the fitting range to am G [0.024,0.030]. Larger quark masses do not seem to fall into the scaling 
region, as the data does not follow the approximate linear behavior for am > 0.032. 

In Fig. H| (right) we show the pseudo critical couplings for different volumes as function of 
mass. For better visibility we have subtracted the linear function j3^ gh (am) = 1.8281- am + 5 .0889 
and multiplied the differences by 10 3 . We see that for large masses (am > am c ) the pseudo crit- 
ical couplings seem to be volume independent within errors, whereas at sub-critical masses we 
observe a distinct scaling with the volume. Assuming that the two masses am = 0.024 and 0.026 
lie in the first order region we perform to extrapolations of the pseudo critical couplings to the 
thermodynamic limit, using f5 pc (N a ) = j6 f (°°) +cN a 3 . Here c is a constant that we fit. We find 
that the infinite volume can, to a good approximation, already be described by the 32 3 -lattices. For 
am < am c we obtain a parameterization of the critical couplings in the thermodynamic limit as a 
function of the quark mass, as /3 c (ara) = 2.12(3) x (am — 0.0268)5.1372(1). This is by definition 
a parameterization of the (negative) t-axis. We can thus obtain a value of the mixing parameter B 
through the relation 

ZH = -^) , (3 .2) 

d ( am ) am=am c 

which arises due to the fact the in the Ising model the line of first order phase transitions defines 
the temperature direction. We estimate B=-0.47(l). Assuming that the scaling directions are or- 
thogonal, this also determines the parameter A. 



4. The mixing parameters r and s 

So far we have determined the critical point as well as the proper scaling directions. Before 
we can consider the magnetic equation of state we still have to clarify the correct order parameter 
M through the determination of the mixing parameter s. We fix s by demanding that the linear 
combination M = yy+sSc fulfills basic properties of an order parameter, namely that M vanishes 
at the critical point and stays zero for t > 0,h = 0. While the former condition is easy to fulfill and 
yields a mixing parameter s = 0.048(2), the latter condition is equivalent with the fact that energy 
and magnetization operators are statistically independent: 

((8M)(8E)) =0 . (4.1) 



Using this relation together with the definition of the scaling operators (2.4) and directions (2.5), 
we obtain the following conditions for the mixing parameters: 

((8 ¥¥ )(SS G ))-B((S ¥¥ ) 2 ) (A0 , 

r = —B and s = —-^ — rr; rm — —rr^ — zr~ ■ (4-2) 

((8S G ) 2 )-B((8 ¥¥ )(8S G )) 

As one can see, we require as input the already determined mixing parameter B, as well as the 
susceptibilities and mixed susceptibilities of the QCD operators ¥¥ and S G . As susceptibilities 
are usually much more difficult do determine, the error on the mixing parameter s is much larger 



using this condition. Moreover, the value for the mixing parameter we obtain from Eq. (12) is 



s = —0.8(1), which is very different form the previously obtained value. We conclude that with 
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Figure 3: The rescaled order parameter M as function of the scaling variable z. 

our ansatz for M we cannot at the same time obtain a very flat behavior for M(t) at t > and 
demand that M(t) vanishes at t = 0. This is most likely due to the finite volume effects that are 
still present in the 24 3 -lattices close to the critical point. We note, that we can obtain a good fit to 
the magnetic equation of state as described in the next section with both values of s. However, as 
long as M does not vanish at the critical point, one is forced to keep the leading order contribution 
from the regular part of the free energy. This will contribute a constant value to M and lead to the 
desired behavior of the order parameter. 



5. Magnetic equation of state 

Close to the critical point, where all regular contributions to the free energy become negligible, 



it follows from ( |2.2| ) that the singular part of the free energy can be expressed as a function of the 
single scaling parameter z = t/h 1 / 8 (by setting the scale parameter b to b = h~ yh ). For the order 
parameter (magnetization), which is the derivative of the free energy with respect to the symmetry 
breaking field h, one then obtains 

M(t,h)=h^ s f G (z). (5.1) 
This relation is known as the magnetic equation of state. The scaling function fc(z) is unique to 



the universality class. In the following we use the parameterization as derived in Ref. Q10|]- In 
order to obtain the last two missing normalization constants to and ho, we fit our data for the order 
parameter M in the range of 0.026 < am < 0.028 to the magnetic equation of state. Using s = 0.048, 
am c = 0.0268 and j3 c = 5.1372 we obtain t Q = 0.0210(9) and h = 0.0005(1). In Fig. | we plot 
the rescaled order parameter M/h 1 ^ 8 as a function of the scaling parameter z. We find good scaling 
behavior within the range of— 10<z<10. 
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6. Outlook 

In our future work we will address the question whether the good scaling behavior of our 
constructed order parameter can also be observed in corresponding susceptibilities. Of special 
interest will be not only the chiral susceptibility but also the mixed susceptibility that contains 
the response to the variation of a small chemical potential (d 2 lnZ/djj,dh = dM/dji). From this 
susceptibility we will be able to determine the curvature of the critical surface that is defined by the 
extension of the critical line in the quark mass plane towards non-zero chemical potentials. 
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